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Abstract 



Q^' The contribution of the sub-process 'jq — s- q'lil2 in hadron-hadron interactions is 



considered. It is a part of one- loop electroweak radiative corrections for the Dreh-Yan 
production of lepton pairs at hadron cohiders. It is shown that this contribution should 
be taken into account aiming at the 1% accuracy of the Drell-Yan process theoretical 
^ ■ description. Both the neutral and charged current cases are evaluated. Numerical results 

m : are presented for typical conditions of LHC experiments. 

PACS: 12.15.Lk Electroweak radiative corrections; 13.40.Ks Electromagnetic correc- 
' tions to strong- and weak-interaction processes; 13.85.Qk Inclusive production with iden- 

. tified leptons, photons, or other nonhadronic particles 

o 

'. 1 Introduction 

> ■ 

><| ■ The Drell-Yan like processes at high energy hadron colliders provide an advanced tool for 

. precision studies of several problems in the elementary particle phenomenology. Studies of 
single Z and W bosons production with the subsequent decays into leptonic pairs play a very 
important role in the physical programs of Tevatron [1,2] and LHC [3,4]. These processes have 
large cross sections and clean signatures in the detectors. That allows to reach at LHC the 1% 
experimental accuracy for the total cross sections of these processes as well as high precision 
in the measurements of differential distributions. In particular, Drell-Yan like processes are 
planned be used at LHC for luminosity monitoring, W mass and width measurement, detector 
calibration, extraction of parton density functions, new physics searches, and other purposes. 

Adequately precise theoretical predictions for single Z and W production at LHC are 
required. For this reason we have to scrutinize several effects involved in the derivation of the 
theoretical accuracy: QCD and electroweak radiative corrections, uncertainties in the partonic 
density functions (PDF's), technical precision of Monte Carlo event generators etc. In this 
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paper we consider a particular contribution of the first order electroweak radiative corrections 
coming from the photon induced process 

hi + h2 X +-f + q ^ X + q' + h + k, (1) 

where hi 2 stand for the initial colliding hadrons; li and I2 is a pair of leptons (e.g. jj~ 
and /i"*", or z/g and e+); X + q' denotes the remaining final state particles (typically they 
are hadrons). Here 7 and q are treated as partons found in the initial hadrons with certain 
energy fractions at a given factorization scale. In this paper we use the MRST2004QED [5] 
parameterization of parton density functions (PDFs), which provides in particular the photon 
content in proton at NLO. Note that the evolution [6] of the partonic densities taking into 
account simultaneous QCD and QED effects leads to the unique value of the factorization 
scale, so that it is impossible to disentangle QED and QCD contributions. This leads also to 
the fact that the reduction of the factorization scale dependence can be reached now only by 
taking into account both QED and QCD higher order radiative corrections and that should 
be performed within the same factorization scheme. Nevertheless due to the smallness of the 
fine structure constant a in comparison with the strong coupling constant a<j, we can limit 
ourselves to the evaluation of only the first order electroweak corrections [7-16] together with 
certain higher order leading logarithmic contributions [17-19]. At the same moment QCD 
corrections have to be treated at least at NNLO [20-22]. Some numerical results for the 
inverse bremsstrahlung contribution to the charged current case (single W boson production) 
were already presented by S. Dittmaier and M. Kramer in the proceedings of the Les Houches 
workshop [23]. We performed an independent calculation and give below a comparison with 
the earlier results. The neutral current case is considered in addition. 

This paper is organized as follows. In the next section we present the derivation of the 
Drell-Yan process cross sections in the scheme with massive quarks. The subtraction of the 
quark mass singularities is described in Sect. [3l Numerical results and their discussion are 
presented in Conclusions. 

2 Inverse bremsstrahlung with massive quarks 

Let us compute the cross section of the process ([T]) in the form proposed by Drell and 
Yan [24] as of a convolution of the parton density functions with the hard sub-process distri- 
bution. In our case the sub-processes is 

q + 1 ^ q' + k + I2, (2) 

where quarks q and q' are of the same type for the neutral current (NC) case and different 
for the charged current (CC) one. We compute the matrix element of the NC and CC sub- 
processes with help of the SANC system [25,26] environment. In the actual calculation, we 
start within the massive quark scheme. The matching of this scheme with the PDF formalism 
will be performed by means of the subsequent subtraction of the quark mass singularities 
from the computed cross section. So, we evaluate the complete tree-level matrix elements of 
the sub-processes in the standard way keeping the exact dependence on the quark and lepton 
masses. The Feynman diagrams for the sub-processes under consideration are shown in Figs. [1] 
andd 



2 




Figure 1: Feynman diagrams for inverse bremsstrahlung in the charged current Drell-Yan 
sub-process. 



We construct the squares of the matrix elements in the usual way and obtain the partonic 
cross sections of the sub-processes. These quantities have to be convoluted then with the 
parton density functions: 

a^inv.brem. ) ^ J dx.dx^Qiix,, M')^{x2, M^ f ^ je(c^, X,, X^) , (3) 

* 

where Ci denotes the cosine of the scattering angle of the first lepton (another variable can be 
chosen as well). The step function 0(ci, Xi, X2) defines the phase space domain corresponding 
to the given event selection procedure. The partonic cross section is taken in the center-of- 
mass reference frame of the initial partons, where the cosine of the first lepton scattering angle, 
Ci, is defined. The transformation into the observable variable Ci involves the Jacobian: 

dc 4:XiX2 , . s 

J = ^ = — 5—, a^Xi+X2-c{Xi-X2), 

dc 

2xi 

c = 1 — (1 — c) , s = SX1X2, (4) 

a 

where s is the squared center-of-mass energy of the colliding hadrons. An analogous formula 
can be written for any other choice of a differential distribution as well as for the total cross 
section. 
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Figure 2: Feynman diagrams for inverse bremsstrahlung in the neutral current Drell-Yan 
sub-process. 

In Eq. ([3]) we presented the contribution, when the photon is found in the first of the 
colhding hadrons and the quark is taken from the other one. Of course, there is also the 
contribution when we choose the particles in the other way round, and it is taken into account 
in our numerical simulations. 



3 Subtraction of the quark mass singularities 

Since the calculation of the partonic cross sections was performed keeping finite masses 
of the quarks, the result depends on the values of the masses. For the high energies 
this dependence arises in the form of large logarithms of the type ln(M^/m^) that give a 
considerable numerical effect, while the other mass-dependent contributions suppressed by 
the factor m^/M^ <^ 1 can be omitted (here M is a typical energy scale of the partonic sub- 
process). The large logarithms represent quark mass singularities. They can be treated with 
help of the QED renormalization group approach. But the point is that they have been already 
taken into account in the evolution of partonic density functions. In the MRST2004QED 
distributions [5] that have been done explicitly. But even in any other PDF on the market 
the QED evolution is implicitly taken into account since it has not been subtracted from the 
experimental data before the PDF fitting procedure. In fact, QED corrections to the quark 
line in deep inelastic scattering are usually omitted in the data analysis, see Refs. [27,28]. 

The quark mass singularity of the first type arises from the right Feynman diagrams in the 
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upper line of Fig. [T] and of Fig. [2] for CC and NC cases, respectively. The singularity originates 
from the kinematical domain when the virtual quark propagator is close to the mass shell. 
For this situation there is a convolution of distributions of the two sub-processes: a conversion 
of the photon into a pair of quarks and a Drell-Yan partonic process g' + g — >■ hh- In the 
MS scheme the corresponding contribution reads 

11 1 



^ii<^i) = j J (ixi dx2 'j{xi, M'^)qi{x2, M'^) j dx3Dg>^{x3, M,m^ 



J Q{ci,XiX3,X2), (5) 







dci 

where ci, J' and s are calculated according to Eq. (jlj) with the interchange xi xix^. For 
the NC case in the above equation we have q' = q. The structure function Dgi^{x3, M,mqi) 
describes the probability to find quark q' with energy fraction x^ in the photon. For the 
MS scheme at NLO this function reads 

D^{x„ M, m,,) = ^Ql \n^[xl + (1 - x.fl (6) 

where M is the factorization scale, and Qg' is the quark charge. 

In the neutral current case there is one additional source of the quark mass singularities. It 
arises from the two lower Feynman amplitudes in Fig. [21 when the virtual photon propagator 
is near the mass shell. In this case we have the convolution of the distributions of the following 
processes: 27 — > /i/i and q — > 7g. The corresponding contribution is 

11 1 
52(ci) = J2j J dxi dx2 gi(xi,M2)7(x2,M2) J dx^D^gix^, M,mg) 

^'0 
X — J Q{ci,XiX3,X2). (7) 

dci 

The relevant structure function describes the probability to find a photon with a certain energy 
fraction in the quark: 

}^JS x3, M,m, = — Q ^ ^ In — 

271 ^ X3 y. 



2?:^;^(X3, M, rUg) = — Q ^ ^ In - 2 Inxs - 1 , (8) 



According to the renormalization formalism we have now to subtract the contributions 
([5]) and (JTj) from the computed cross section ([3]). In a realistic situation we have to perform 
this procedure numerically in order to keep the possibility to impose experimental cuts. On 
the other hand, it can be shown analytically that the terms with the logarithms of the quark 
masses do cancel out during the subtraction procedure. 
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4 Numerical Results and Conclusions 



For the numerical evaluations we used the same conditions and the input parameters as in 
Ref. [23]: 



G p 


= 1.16637 X 10-^ GeV~^ 






a(0) 


= 1/137.03599911, 




= 0.1187, 




= 80.425 GeV, 


r 

J- w 


= 2.124 GeV, 


M, 


= 91.1867 GeV, 




= 2.4952 GeV, 




= 150 GeV, 


mt 


= 174.17 GeV, 


m„ 


= = 66 MeV, 




= 1.55 GeV, 


ITLs 


= 150 MeV, 


rrib 


= 4.5 GeV, 


\Vud\ 


= |Ks| = 0.975, 


\Vus\ 


= |Krf| = 0.222 



The MRST204QED set [5] of PDF's and the Gp EW scheme were used. The cut on the 
charged lepton rapidity and transverse momentum are jry^l < 1.2 and Pt/ > 25 GeV. The cut 
on the missing transverse momentum for the CC case is imposed as well: -Pr.missing > 25 GeV. 

At the partonic level for the CC and NC processes (7 + g q' + h + h) we performed a com- 
parison with the corresponding distributions obtained with help of the CompHEP system [29] 
and found a good agreement. 





25-00 


50-00 


100 - 00 


200 - 00 


500 - 00 


1000 - 00 


cjo/pb 


DK 


2112.2(1) 


13.152(2) 


0.9452(1) 


0.11511(2) 


0.0054816(3) 


0.00026212(1) 


SANG 


2112.2(1) 


13.151(1) 


0.9451(1) 


0.11511(1) 


0.0054813(1) 


0.00026211(1) 




DK 


0.071(1) 


5.24(1) 


13.10(1) 


16.44(2) 


14.30(1) 


11.89(1) 


SANG 


0.074(1) 


5.24(1) 


13.09(1) 


16.43(1) 


14.30(1) 


11.90(1) 



Table 1: Cross sections ctq and a^q of the processes p[q]p[q'] — > z/^yU"^X and p[7]p[q] z/^/i"'"X, 
respectively and corresponding corrections = cr^q/ao, obtained by DK and SANG groups 
for different Pt,^ ranges at LHC. 

In Tabled] we present the results of comparison for the inverse bremsstrahlung contribution 
to the CC Drell-Yan process with different cuts on the charged lepton transverse momentum 
(see the details in Ref. [23]). Our results are marked as "SANG", they are compared with 
the numbers ("DK") presented by the S. Dittmaier and M. Kramer in Ref. [23]. The small 
deviations in the results for the values of the corrections are certainly beyond the 1% precision 
level. They are due to some differences in the schemes of calculations and are induced by 
higher order effects in a. 

Table [2] shows the results of comparison for the inverse bremsstrahlung contribution to 
CC Drell-Yan process with different cuts on transverse mass of muon-neutrino pair. The 
corresponding numbers for are below percent level. 

Table [3] gives the results for the inverse bremsstrahlung contribution to the neutral current 
Drell-Yan process with production of two muons. Different values of the cut on the invariant 
mass of the muon pair are considered. For the Born cross section we show also the numbers 
of HORACE [17, 19], which are in fair agreement with the SANG results. 
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50-00 


100 - oo 


200 - oo 


500 - oo 


1000 - oo 


2000 - oo 


cjo/pb 


DK 


2112.2(1) 


13.152(2) 


0.9452(1) 


0.057730(5) 


0.0054816(3) 


0.00026212(1) 


SANG 


2112.2(1) 


13.151(1) 


0.9451(1) 


0.057730(5) 


0.0054813(1) 


0.00026211(1) 




DK 


0.0567(3) 


0.1347(1) 


0.2546(1) 


0.3333(1) 


0.3267(1) 


0.3126(1) 


SANG 


0.0532(1) 


0.1350(1) 


0.2537(1) 


0.3314(1) 


0.3245(1) 


0.3094(1) 



Table 2: Cross sections cxo and a^q of the processes p[q]p[q'] I'^^'^X and p[7]p[q] u^^fi'^X, 
respectively and corresponding corrections = cr^q/o"o, obtained by DK and SANC groups 
for different Mt,,/^^+ ranges at LHC 



M^+,,-/GeV 


50-00 


100 - oo 


200 - oo 


500 - oo 


1000 - oo 


2000 - oo 


cro/pb 


HORAGE 


254.64(1) 


10.571(1) 


0.45303(3) 


0.026996(2) 


0.0027130(2) 


0.00015525(1) 


SANG 


254.65(2) 


10.571(1) 


0.45308(3) 


0.026996(2) 


0.0027131(2) 


0.00015525(1) 




SANG 


0.047(1) 


0.449(1) 


0.013(1) 


0.496(1) 


0.619(1) 


0.563(1) 



Table 3: Cross sections do and cr^q of the processes p[q]p[q'] — ^ fi'^fi X and p[7]p[q] — ^ fi'^fJ^ X, 
respectively and corresponding corrections = cr^q/cxo, for different M^+^- ranges at LHC. 



In Fig. [3] we plotted the distributions of the Born-level cross section and of the relative 
radiative correction versus the transverse mass of the muon and neutrino pair MT(/i"'"i^^) in 
the CC Drell-Yan process, 

MT(/i+Z/^) = ^J2PT,^PTM-COS(j)^,), (9) 

where is the angle between the muon momentum and the missing one in the transverse 
plane. In Fig. H] the analogous distributions in the muon transverse momentum P^.^i are given. 

Fig. [5] shows the Born differential cross section of the neutral current Drell-Yan process 
(a) and the relative correction S^q (b) as a function of invariant mass M^+^- of the muon pair. 
Fig. O gives us results for the Born differential cross section of the neutral current Drell-Yan 
process (a) and the relative correction 6^q (b) as a function of /i"*" transverse momentum Pt,/^- 
The distributions around the W and Z resonances are plotted. The drop-offs in the first 
bins of the correction distributions in NC have no any physical sense. They arise because the 
factorization procedure with the longitudinal partonic density functions doesn't allow to apply 
the experimental cuts unambiguously. The drop-offs can be shifted by choosing a different cut 
value. We checked that the rest of the distributions doesn't suffer from this problem. 

In this way we presented the photon-induced contribution to the first order electroweak 
radiative corrections to Drell-Yan processes. For the case of charged current scattering our 
results are in a good agreement with earlier calculations of the other group. The neutral 
current case was considered in an analogous manner. This inverse bremsstrahlung contribution 
should be taken into account together with all other relevant effects to reach the accuracy of 
the Drell-Yan process theoretical description adequate to the precision of the forthcoming 
LHC experiments. The typical size of the contribution is below one percent, but for the 
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Figure 3: The Born-level CC Drell-Yan cross section and the relative contribution of the 
inverse bremsstrahlung versus the transverse mass of the muon-neutrino pair. 




25 30 35 40 45 50 25 30 35 40 45 50 



Pj(p.') [GeV] Pjiti-) [GeV] 

(a) (b) 



Figure 4: The Born-level CC Drell-Yan cross section and the relative contribution of the 
inverse bremsstrahlung versus the /i^ transverse momentum. 

case of transverse momentum distribution in CC scattering, the effect can reach up to 16% 
depending on the cut value. We are going to implement the results of our calculations into a 
general Monte Carlo event generator for Drell-Yan processes, which is under development in 
the SANC group. 
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Figure 5: The Bom- level NC Drell-Yan cross section and the relative contribution of the 
inverse bremsstrahlung versus the invariant mass of the muon pair. 
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Figure 6: The Born- level NC Drell-Yan cross section and the relative contribution of the 
inverse bremsstrahlung versus the /i^ transverse momentum. 
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